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A  Simulation  Analysis  of  Errors  in  the 
Measurement  of  Standard  Electrochemical  Rate  Constants  from 


Phase -Selective  Impedance  Data 


David  F.  Milner  and  Michael  J.  Weaver 
Department  of  Chemistry 
Purdue  University 

West  Lafayette,  Indiana  47907  (U.S.A.) 

ABSTRACT 

A  digital  simulation  analysis  is  presented  of  the  coupled  effects  of 
solution  resistance  and  lock-in  amplifier  damping  upon  the  evaluation  of 

g 

standard  rate  constants,  k^,  f°r  rapid  electrode  reactions  using  phase - 
selective  AC  impedance  measurements  with  positive-feedback  iR  compensation. 
These  two  effects  combine  so  to  yield  "apparent"  (i.e.  measured)  rate 

g 

constants,  kQlj(app),  evaluated  using  the  conventional  frequency -dependent 
analysis  that  are  generally  smaller  than  the  actual  values,  k^(true).  The 
extent  to  which  these  systematic  errors  depend  on  experimental  conditions 

g 

is  explored  for  various  values  of  ko^(true),  the  uncompensated  and  specific 
solution  resistances,  the  amplifier  filter  time  constant,  the  double -layer 
capacitance,  and  other  relevant  parameters.  A  simple  scheme  employing 

g 

these  simulations  is  outlined  by  which  k^Ctrue)  values  can  be  extracted 

from  kj^app)  measurements.  Some  experimental  verification  of  these 

effects  is  also  presented.  Accession  For 
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The  measurement  of  rate  constants  for  electron  exchange  of  redox 
couples  at  metal -solution  interfaces  (so-called  "standard"  electrochemical 
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rate  constants,  k^)  has  long  been  a  matter  of  fundamental  interest.  A 
large  and  important  class  of  reactions  involve  redox  couples  for  which  the 
inner-shell  barrier  (i.e.,  that  arising  from  structural  differences  between 
the  redox  forms)  is  negligible  or  small,  so  that  the  activation  energy  AG* 
arises  predominantly  from  outer-shell  (solvent)  reorganization.  Current 
interest  in  such  processes,  as  well  as  in  self-exchange  reactions  in 
homogeneous  solution,  centers  on  the  anticipated  important  contribution  of 
solvent  dynamics  to  the  preexponential  factor.  For  such  systems  the  k  ^ 
is  predicted  to  approach  and  possibly  surpass  the  maximum  values  amenable 
to  experimental  evaluation,  even  when  using  rapid  perturbation  techniques. 
Consequently,  under  these  circumstances,  estimates  of  obtained  using 
such  methods  may  be  entirely  erroneous.  It  is  therefore  of  paramount 
importance  to  both  identify  and  quantify  the  various  experimental  factors 
that  determine  this  upper  limit  when  using  a  given  instrument  and 
electrochemical  technique . 

A  primary  method  used  to  evaluate  such  rapid  electrochemical  rate 
constants  involves  AC  impedance  measurements  with  phase -selective 
detection. ^ The  contributions  of  Don  Smith  in  this  area  are  well  known; ^ 
his  work  has  had  a  central  influence  on  the  development  of  such  AC  methods 
and  their  broadbased  application  in  electrochemistry,  including 
applications  to  electrode  kinetics.  The  cell  impedance  as  a  function  of 
the  AC  frequency  and  the  electrode  potential  can  provide  a  particularly 


straightforward  route  to  the  determination  of  k^-  The  technique  as  it 
is  commonly  practiced  utilizes  a  potentiostat  with  positive-feedback  IR 


VC* 
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lock-in  amplifier  employed  as  a  phase -sensitive  detector,  and  either  a 
dropping  mercury  electrode  (DHE)  or  a  stationary  mercury  or  solid  surface 
for  the  working  electrode. 

Of  the  various  factors  that  can  conspire  to  produce  a  significant  and 

5 

ultimately  limiting  impediment  to  the  reliable  determination  of  k^,  the 
influence  of  solution  resistance,  Rs,  is  paramount.  Besides  the  common  use 
of  positive-feedback  iR  compensation,  the  presence  of  Rs  is  accounted  for 
in  a.c.  impedance  measurements  by  means  of  a  vectorial  subtraction 
scheme,  originally  due  to  Randles.^  However,  in  both  these  approaches  the 

g 

very  presence  of  R_  will  place  inevitable  limitations  on  the  values  of  k  , 

ob 

that  can  be  reliably  evaluated,  especially  for  large  Rs  values.  Analyses 

dealing  with  the  manner  and  extent  to  which  these  limitations  occur  in 

practice  have,  nevertheless,  remained  surprisingly  sparse. 

We  present  here  in  some  detail  a  digital  simulation  analysis  of  this 

problem,  primarily  as  applied  to  a.c.  impedance  measurements  employing 

positive -feedback  iR  compensation.  We  have  briefly  reported  previously  a 

simple  approach  which  can  enable  estimates  of  the  error  induced  in  kS, 

ob 

measurements  using  the  conventional  analysis,  based  on  simulations  using 
trial  values  of  Rs  and  the  double-layer  capacitance,  The  present 

report  contains  a  development  of  this  analysis,  and  also  provides  estimates 
of  the  additional  measurement  error  that  can  be  induced  by  the  solution 
resistance  when  coupled  with  damping  in  the  lock-in  amplifier.  Some 
experimental  verification  of  these  effects  is  also  presented.  The  origin 
and  magnitude  of  the  errors  from  measuring  the  AC  currents  themselves  that 
result  from  amplifier  damping  are  detailed  in  a  separate  paper. ^ 

The  overall  aim  of  this  work  is  to  aid  the  experimentalist  in 
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ascertaining  the  validity  as  well  as  the  quantitative  accuracy  of  fast 
electrode  kinetic  measurements  made  using  phase -selective  a.c.  detection. 
Taken  together,  we  believe  that  these  analyses  provide  a  relatively 
straightforward  means  by  which  the  influence  of  these  undesired  solution 
resistance  and  amplifier  damping  effects  can  be  diagnosed  as  well  as 
minimized. 

Simulation  Procedures 

The  general  approach  employed  here  involves  producing  digital 
simulations  of  the  AC  current,  including  the  time  dependence  at  a  growing 
DME,  at  a  given  fixed  potential  either  in  the  presence  or  the  absence  of  an 
electroactive  species. 

We  employ  the  conventional  equivalent  circuit  model  of  an 
electrochemical  cell  as  it  is  formulated  for  simple  redox  processes  in  the 
AC  experiment,  shown  schematically  in  Fig.  1.  (This  model,  originally  due 
to  Randles,5  is  generally  accepted  as  a  reasonable  approximation.8)  It  is 
readily  apparent  that  the  cell  impedance  is  given  by: 

zcell  “  <  1  +  1  )  +  Rus  U> 

ZCdl  RF  +  ZCF 

where  Rug  is  the  uncompensated  solution  resistance  (i.e.  that  part  of  the 
total  system  resistance  which  has  not  been  compensated  for  electronically) , 
Rp  is  the  resistive  part  of  the  faradaic  impedance,  Zq^  is  the  impedance 
of  the  double-layer  capacitance,  and  Zqf  is  the  impedance  of  the  capacitive 
part  of  the  faradaic  impedance. 

The  Impedance  of  a  capacitance  is  given  by: 


where  w  is  the  radial  frequency,  and  j  is  the  unit  irrational  constant 

(T7!)- 


For  the  DME,  one  must  consider  the  time  dependence  of  the  electrode 
area,  A.  Besides  the  obvious  proportionality  between  and  A,  the 

solution  resistance,  and  thus  the  uncompensated  resistance,  of  an 
electrochemical  cell  is  also  a  function  of  A  provided  that  the  level  of 
compensation  is  not  altered  as  the  area  varies.  For  a  spherical  electrode, 
such  as  is  approximated  by  a  DME,  the  uncompensated  resistance  as  a 
function  of  surface  area  is  given  by:^ 

®us  “  ^hom  +  P (^*A)  **  (3) 

where  p  is  the  specific  solution  resistance,  Rhom  *s  t*'e  "homogeneous"  part 
of  the  resistance  (that  portion  which  is  independent  of  the  electrode  area) 
and  Rus  is  the  uncompensated  solution  resistance.  Since  the  term  p(4w A)*15 
represents  the  resistance  between  a  spherical  electrode  and  a  reference 
electrode  in  the  bulk  solution,  Rhoin  is  approximately  zero  in  the  absence 
of  iR  compensation. ^ 

Equations  describing  Rp  and  Cp  resulting  from  the  presence  of  the 
electroactive  species  are  also  necessary.  These  can  be  expressed  as:^a 

Rf  -  (wCp)-1^  +  (2<A/A)]  (4) 

where 

X  -  (k*  /D**)(exPr-  °nF  (E-E*)]  +  exp  {(1-a)  (E-E°)]j  (4a) 

RT  .  RT 

r  wn2F2AC*(2wD)li 

°F  “  - - -  (5) 

4RT  coshz(RT/2nF)(E-E8) 

In  these  equations,  D  is  the  reactant  diffusion  coefficient,  a  is  the 
cathodic  transfer  coefficient,  n  is  the  number  of  electrons,  E  is  the 
electrode  potential,  E°  is  the  standard  potential  of  the  redox  couple,  and 


5 


A  is  the  electrode  area. 

In  the  DME  case,  the  time  dependence  of  A  is  given  by: 
A  -  4w(3mt/4iro)2  /  3 


(6) 


where  m  and  a  are  the  mercury  flow  rate  and  density,  and  t  is  the  time 
since  drop  birth. 

These  relations  enable  us  to  digitally  simulate  in-phase  and 
quadrature  currents  as  a  function  of  (E-E*)  and  the  applied  frequency  w, 
for  various  input  parameters,  notably  k^.  Rs  (&uS  when  iR  compensation  is 
involved),  and  C^.  For  simplicity,  parameters  which  have  a  minor  or 
easily  predictable  impact  upon  the  results  were  held  constant  at  the 
following  typical  values:  D  -  10"^  cm^s"^,  C*  -  1  mj;  ,  m  -  3  mg  s'*,  tD 
(the  mercury  drop  time)  -  2s,  T  -  300K,  n  -  1,  a  -  0.5,  and  EAC  (the  rms 
amplitude  of  the  applied  AC  signal)  -  2.24mV. 

These  simulated  impedance  data  were  then  used  to  obtain  "apparent" 
values  of  k®^>  k^(app) ,  using  the  conventional  analysis  involving  the 
frequency -dependence  of  the  phase  angle.  -*a  Specifically,  the  cotangent  of 
the  phase  angle,  cot  was  taken  to  be: 


.  ,  _cor/Tcor 

cot  4>  -  IQ  /Ij 


(7) 


where  Ip  and  Ij  are  the  simulated  quadrature  and  in-phase  currents 
corrected  for  the  corrected  background  currents  sufficiently  far  from  the 
AC  polarographic  wave  so  that  faradaic  currents  are  negligible.  The  cot  ^ 
values  were  then  plotted  against  w1*,  and  the  slope  of  the  resulting 
regression  line,  s,  used  to  obtain  k^(app)  from^a: 


where 


i 

G  -  (exp  a)/(l  +  exp  a)  (8a)  j 

g 

The  quantity  kQ^(app)  is  therefore  the  rate  constant  that  would  be 
extracted  from  AC  impedance  measurements  with  an  electrochemical  cell  , 

t 

having  equivalent  circuit  components  equal  to  the  simulation  input  ' 

parameters.  The  extent  to  which  k^(app)  differs  from  the  "true"  value, 

< 

s 

kQ^(true) ,  represents  the  systematic  error  which  will  necessarily  occur 

when  the  conventional  analysis  is  employed  under  a  given  set  of 

experimental  conditions.^  In  terms  of  the  present  analysis,  two  factors 
s  s 

will  cause  kQ^(app)  to  differ  from  kQ^(true).  Firstly,  when  the 
uncompensated  solution  resistance,  Rug,  is  not  exactly  zero,  the 
conventional  analysis  embodied  in  Eq.  (7) -(8)  is  not  strictly  valid.  The 
extent  of  the  resulting  error  clearly  depends  upon  Ryg ;  the  simulation 
results  given  below  provide  a  means  of  estimating  the  magnitude  of  the 
correction  to  k^(app)  that  is  required. 

g 

An  additional  source  of  error  in  determining  k  ^  if  significant  time- 
dependent  currents  are  involved  (most  commonly  when  using  a  DME)  arises 
from  the  distortion  in  the  measurement  of  time -dependent  in-phase  and 
quadrature  currents  using  a  lock-in  amplifier  that  results  from  a  damping 
effect.  The  distortion  of  the  signal  caused  by  the  bandwidth- limiting 
portion  of  the  lock-in  amplifier,  namely  the.  low-pass  filter  stages 
following  the  balanced  mixer,  can  be  described  by^>^® 

t 

V(t)  -  (RC)’1/  Vi(r)exp[-(t-r)/RT]dr  (9) 


where  V(t)  and  Vi(t)  are  the  voltage  at  the  filter  output  and  input, 


respectively,  at  time  t,  and  RC  is  the  filter  time  constant.  It  is 
apparent  from  Eq.  (9)  that  for  the  simple  RC  low-pass  filters  used  in  most 
commercial  lock-in  amplifiers,  the  resulting  absolute  error  in  the  measured 
current  is  approximately  proportional  to  the  time  derivative  of  the 
detected  signal.  Further  details  of  this  aspect  of  the  simulations  are 
given  elsewhere. ^  Such  damping  effects  turn  out  to  be  of  significance  even 
when  Rus  -  0  if  relatively  high-resistance  media  (i.e.  high  p)  are 
employed. 

RESULTS  AND  DISCUSSION 

The  uncompensated  solution  resistance,  R^,  may  be  regarded  as  being  a 
resistive  component  in  series  with  the  electrochemical  double  layer.  If 
the  AC  cell  current  is  non-zero,  this  resistance  will  in  general  prevent 
the  phase  of  the  double -layer  potential  from  equalling  that  of  the  AC 
potential  applied  to  the  cell.  If  Rus  is  positive,  the  phase  angle  of  the 
double-layer  potential  will  be  between  0’  and  -90°  with  respect  to  the 
applied  potential  In  AC  voltammetric  measurements.  Since  the  phase  angle 
of  the  double-layer  potential  is  negative,  the  measured  phase  angle  of  the 
faradaic  current  will  then  be  less  than  that  with  respect  to  the  actual 
double -layer  potential.  This  results  in  larger  values  for  cot  a 
greater  regression  slope,  and  thus  smaller  values  of  k®b(app)  than 
k®b(true).  We  will  now  examine  quantitatively  the  effects  of  positive  R^ 
values,  as  these  are  predicted  to,^  and  in  our  laboratory  do,  arise  more 
frequently  than  negative  values  for  Rus. 

Effects  of  Uncompensated  Resistance  at  a  Stationary  Electrode 


The  case  of  a  stationary  electrode  will  be  considered  first  since  the 


effects  of  lock-in  amplifier  damping  will  usually  then  be  absent.  In 
addition  to  the  fixed  simulation  parameters  noted  above,  we  set  A  -  0.028 


8 


cm^,  which  is  the  area  obtained  for  a  DME  at  the  end  of  drop  life  using  the 
previously  quoted  values  for  m  and  tj>.  Simulations  are  presented  here  for 
two  (typically  "small"  and  "moderate")  values  of  the  double-layer 
capacitance,  C^,  5  and  20  pT  cm*^.  In  the  absence  of  amplifier  damping 
effects,  the  magnitude  of  the  total  solution  resistance  will  not  directly 

affect  the  simulation  results,  since  only  the  uncompensated  component 

s  s 

will  cause  kQ^(app)  to  differ  from  kQ^(true). 

It  is  useful  to  consider  simulations  aimed  at  estimating  the  upper 

limit  to  k^  that  can  be  obtained  for  a  particular  electrode  -  electrolyte 

arrangement,  as  dictated  by  the  minimum  R^  value  that  can  be  achieved. 

This  is  most  simply  gleaned  from  the  values  of  k®b(app),  which  we  label 

k®b(app,lim) ,  that  are  extracted  from  simulated  data  with  k®b(true) 

taken  as  infinity  (i.e.,  immeasurably  large).  The  experimental  observation 
s  s 

of  kQb(app)  values  that  approach  kob<app,lim)  for  a  particular 

experimental  arrangement  provide  a  clear  diagnosis  of  significant  errors 

in  kQb(app)  arising  from  a  non-negligible  Kus  value. 

Curve  A  of  Fig.  2  is  a  plot  of  k^b(app,lim)  versus  log  R^g  in  the 

absence  of  amplifier  damping  effects,  and  for  Cjji  -  5ftF  cm*^.  This  plot 

clearly  shows  that  even  a  small  value  of  uncompensated  solution  resistance 

can  be  a  serious  impediment  to  the  measurement  of  kS,  for  a  rapid  reaction. 

ob 

For  example,  the  k®b(app,lim)  value  with  only  1  n  of  uncompensated  solution 
resistance  is  15  cm  s*^-.  For  a  more  typical  value  of  3  ohms  for  R^ ,  such 
as  is  commonly  attainable  using  a  PAR  173/179  potentiostat ,  for  example, 
the  value  of  k*b  (app.lim)  decreases  to  5  cm  s'*.  Even  larger  values  of 


'Rvls  are  to  be  expected  experimentally  as  one  moves  to  less  polar,  and  hence 

more  resistive  solvents,  with  the  increased  possibility  that  the  intended 
s 

measurement  of  k  ^  will  contain  information  only  on  the  magnitude  of  R^ 
rather  than  on  the  reaction  rate. 

Curve  A  of  Fig.  3  illustrates  the  corresponding  relationship  between 
k  ,  (app.lim)  and  log  for  -  20  fiF  cm"^ .  While  this  curve  appears 
to  have  the  same  qualitative  shape  as  curve  A  of  Fig.  2,  the  decrease  in 

g 

k  ,  (app.lim)  with  increasing  is  significantly  steeper  in  Fig.  3.  Over 

g 

the  R^g  range  examined,  0.5  to  100  fl,  the  values  of  k^^app.lim)  for 
Cdi  -  5  /iF  cm" ^  are  about  twice  the  corresponding  value  for  -  20  fiF 

cm'2.  This  increase  in  the  degree  of  distortion  for  larger  is  due 

solely  to  the  increased  nonfaradaic  AC  current  which  must  flow  through  the 
uncompensated  solution  resistance,  causing  a  greater  phase  shift  in  the 
double -layer  potential. 

In  order  to  provide  a  more  quantitative  estimate  of  these  errors,  it 

s  s 

is  useful  to  examine  the  relation  between  kQ^(app)  and  kQ^(true)  for 

particular  Rus  values.  Curve  A  of  Fig.  4  is  a  plot  of  log  k^(app) 
s 

versus  kQ^(true)  for  Rus  -  2  11,  again  in  the  absence  of  amplifier  damping. 

g 

As  might  be  expected,  it  is  seen  that  the  relative  error  in  kQ^(app)  is 
relatively  small  for  k^(true)  -  0.5  cm  s"^-,  and  increases  rapidly  for 

g 

larger  rate  constants  as  kQ^  reaches  the  asymptotically  limiting  value, 
kob(app,llm) ’ 

Having  exposed  the  degree  to  which  relatively  small  R^  values  can 

g 

cause  substantial  errors  in  k^  measurements  using  the  conventional 
analysis  based  on  the  dependence  of  cot^  upon  w1*,  it  is  reasonable  to 
inquire  if  such  errors  can  be  diagnosed  from  the  cot^-w^  plots  themselves. 


That  this  is  indeed  the  case  is  illustrated  in  Fig.  5,  which  shows 


simulated  cot^-w1*  plots  obtained  for  several  values  of  R^g  (0,  0.5,  1, 
and  2  0),  again  with  kQ^(true)  -  <®.  The  salient  feature  of  Fig.  5  is  that 
progressively  more  positive  values  yield  not  only  increasing  cot^-w*1 
slopes,  but  also  produce  increasingly  sub-unity  y-intercepts  and  a 
noticeable  upward  curvature  in  the  plots.  These  features  can,  at  least  in 
principle,  be  used  to  distinguish  between  nonzero  slopes  arising 

from  nonzero  R^g  values  and  from  finite  electrode  kinetics  since  the  plots 
should  always  be  linear  with  a  y-intercept  of  unity  if  the  latter,  desired, 
feature  dominates . 

Solution  Resistance  Effects  at  a  DME 

As  noted  above,  at  a  non- stationary  electrode  (most  commonly  a  DME) 

where  the  time  derivative  of  the  amplitude  of  the  AC  current  is  non-zero, 

the  measurement  error  can  be  greater  than  It  would  be  at  a  stationary 

electrode  under  similar  conditions  due  to  the  necessary  low-pass  filter 

present  on  the  output  of  the  lock-in  amplifier. ^ 

Figures  2  and  3  also  contain  a  series  of  plots  of  log  k^(app.lim) 

versus  log  Ryg  for  progressively  increasing  levels  of  low-pass  filtering. 

These  were  computed  for  a  specific  solution  resistance,  p,  of  10  Q  cm, 

which  corresponds  here  to  Rs  -  170.  (This  p  value  is  characteristic  of 

concentrated  aqueous  electrolytes.)  Besides  curve  A  in  these  figures, 

which  as  noted  above  refers  to  no  filtering,  curves  B,  C,  and  D  refer  to 

0.03  s,  0.1  s  time  constant,  and  two  0.1  s  time  constants  in  series, 

respectively.  They  show  that  the  additional  error  in  the  measurement 

induced  by  low-pass  filtering  is  essentially  negligible  for  Rus  ^20. 
s 

As  ko^(true)  is  decreased  to  measurable  values  this  error  becomes  even 


> 


smaller.  Generally  speaking,  then,  we  can  conclude  that  the  additional 
error  due  to  reasonable  levels  of  low-pass  filtering  can  normally  be  safely 
neglected  in  such  highly  conducting  electrolytes. 

For  less  conducting  solvents  such  as  most  nonaqueous  media,  however, 
the  situation  is  quite  different.  Figure  6  consists  of  plots  of 

g 

log  kQ^(app)  versus  log  R^  computed  for  a  markedly  higher  specific 

resistance,  p  -  500  fl  cm;  other  chosen  parameters  are  kQ^(true)  -  1  cm  s'1 

and  C^i  -  20  pF  cm*^.  For  R^  £  2.5  0,  the  error  in  k^(app)  is  less  than 

25%  in  the  absence  of  low-pass  filtering.  However,  in  the  presence  of  as 

little  distortion  as  provided  by  a  single  filter  with  a  30  ms  time 

if  Ryg  -  0  this  error  becomes  nearly  40%,  and  when  R^  -2.50  the  value  of 
s  s 

kQ^(app)  becomes  only  one  half  of  kQ^<true) .  Since  the  extent  of  the  error 
caused  by  the  lock-in  damping  is  approximately  proportional  to  the 

magnitude  of  the  currents  involved,  these  errors  will  decrease  for  smaller 

s  s 

values  of  kQb(true).  For  example,  with  an  Rus  of  1  0,  a  kQb(true)  of 

1  s 

0.5  cm  s'A  and  a  30  ms  filter,  the  error  in  kQb(app)  would  be  only  30%. 

These  effects  of  lock-in  amplifier  damping  are  also  seen  in  the  log 
s  s 

log  kQb(app)  -  log  kQb(true)  plots  given  in  Fig.  4;  curves  B  to  D  refer  to 
the  presence  of  increasing  damping,  for  the  same  conditions  as  in  Fig.  6. 

Decreasing  the  double-layer  capacitance  to  5  pF  cm*^,  while 
maintaining  the  specific  resistance  at  500  n  cm,  produces  significant 
changes  in  the  effects  of  both  lock-in  damping  and  uncompensated 
solution  resistance.  In  the  absence  of  low-pass  filtering,  for  k^(true)  - 
1  cii  s'1,  up  to  5  0  of  Rys  may  be  present  before  the  error  in  kQb(app) 
reaches  25%.  Vith  a  DME,  on  the  other  hand,  for  a  30  ms  time  constant 
filter  is  used  only  1  fi  of  R^  needs  to  be  present  before  a  25%  error  in 


k^(app)  is  encountered.  Again,  as  ko^(true)  decreases  the  situation 
becomes  progressively  improved,  so  that  for  kQ^(true)  -  0.5  cm  s'-1-,  5  0  or 
less  of  uncompensated  solution  resistance  results  in  less  than  a  25%  error 

g 

in  k  ,  (app),  even  when  a  30  ms  filter  is  employed. 


The  foregoing  simulations  demonstrate  that  it  is  desirable  to 
minimize,  or  at  least  have  an  accurate  knowledge  of,  in  order  to 
evaluate  rapid  rate  constants  with  the  conventional  analysis  using  AC 
impedance  measurements.  In  principle,  the  simplest  method  of  estimating 
Ryg  is  from  the  cell  impedance  either  at  potentials  distant  from  the  AC 
polarographic  wave,  or  in  the  absence  of  the  electroactive  species  but 
using  identical  iR  compensation.  In  practice,  however,  the  measurement  of 
Rus  to  the  desired  accuracy  (ca.  1-2  Q)  is  difficult  if  not  impossible  in 
this  fashion.  For  example,  if  Cd^  -  20  pF  cm"^  we  require  there  to  be  3  [i 
of  uncompensated  resistance  before  Ij  becomes  1%  as  large  as  Iq  even  at 
1000  Hz.  Operating  at  lower  frequencies  or  with  smaller  Cd^  values 
decreases  this  fraction  even  further.  Extremely  precise  Ij  measurements 
are  clearly  needed  even  to  approximately  measure  R^.  Thus  under  these 
conditions,  about  a  2%  change  in  Ij  (equivalent  to  one  degree  of  phase 
shift)  corresponds  to  about  5  0.  Such  measurements  are  further  complicated 
by  the  presence  of  significant  phase  distortion  in  most  potentiostats ;  this 
places  an  effective  upper  limit  on  the  frequencies  that  can  be  employed. 
(We  have  examined  a  PAR  174/50,  a  PAR  173/179,  a  PAR  273,  and  a  Hi-Tek 
DT2101  potentiostat  and  found  them  all  to  suffer  somewhat  in  this  regard.) 

We  have  described  one  approach  which  at  least  partly  circumvents  this 
difficulty.^  This  involves  employing  an  electronic  dummy  cell  with 


components  (C^,  Rs)  that  approximately  match  those  of  the  actual  cell 

being  employed.  After  achieving  the  optimal  level  of  iR  compensation,  the 

Cdi  value  is  increased  sufficiently  (5  to  20  fold)  so  that  readily 

measurable  Ij  values  could  be  obtained,  thereby  yielding  Rus  from  the  usual 

2  2 

formula  R^  -  Ea^,[Ij/(Iq  +  1^)].  (The  necessary  assumption  is,  of  course, 
that  Rus  is  independent  of  . ) 

A  number  of  authors  have  evaluated  rate  constants  for  fast  electrode 
reactions  with  AC  impedance  in  the  presence  of  large,  ;>  100  0,  values  of 
the  uncompensated  resistance  by  measuring  the  value  of  R^g  and  subtracting 
it  from  the  cell  response  by  using  the  vectorial  scheme  first  described  by 
Randles.-*  This  approach,  however,  apparently  provides  no  benefits  over  the 
use  of  maximal  iR  compensation  as  considered  here .  Thus  even  though  the 
measurement  of  large  values  of  Rus  can  be  made  with  a  greater  relative 
precision  (AR/R)  than  for  small  values,  the  absolute  precision  is  no 
better.  It  is  clearly  the  latter  which  is  important  in  calculating  the 
error  in  k  ,  (app)  arising  from  the  solution  resistance. 


As  in  practice  it  is  nearly  impossible  to  escape  the  presence  of  at 
least  a  small  amount  of  uncompensated  solution  resistance,  the  effects  of 
lock-in  damping  upon  measured  rate  constants  cannot  be  examined  experi¬ 
mentally  in  the  absence  of  the  influence  of  Rus-  Also,  as  the  most 
reliable  rate  measurements  are  made  using  a  DME  where  problems  of  surface 
contamination  are  minimized,  we  chose  to  examine  both  effects  together. 

Our  approach  was  to  select  a  redox  couple,  cobalticinium-cobaltocene 
(CP2C0"4"/0) ,  that  exhibits  relatively  well-defined,  yet  measurable, 
electrode  kinetics  at  the  DME  in  a  number  of  nonaqueous  solvents. ^  By 
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using  low  mercury  flow  rates,  long  drop  times,  minimal  lock-in  filtering, 


and  maximal  iR  compensation  the  corrections  that  are  necessary  to  apply  to 


k8fe(app)  so  to  yield  k*b(true)  can  be  arranged  to  be  moderate  or  minor. 


For  example,  in  dime thylsulf oxide  (DMSO)  containing  0,1  U  tetrabutyl- 
ammonium  hexafluorophosphate  we  obtain  kflb(true)  -  0.7  cm  s"  ,  the  corres¬ 


ponding  ks,  (app)  values  being  only  20-30%  smaller. 


Using  a  DME  of  higher  flow  rate,  2.6  mg  s‘^,  and  a  larger  amplifier 


time  constant,  0.1  s,  we  then  determined  k  ,  (app)  as  a  function  of  Rug. 


To  do  this,  we  obtained  values  of  kS.  (app)  for  Cp2Co+/°  at  maximal  iR 


compensation  (i.e.  minimal  R^,  at  the  setting  just  prior  to  potentiostat 


oscillation)  and  then  diminished  the  compensation  level  to  be  10,  20,  and 


then  50  O  below  this  point.  (See  ref.  12  for  experimental  details.)  A 


comparison  between  the  resulting  k  ,  (app)  values  and  those  simulated  from 


the  additional  Rus  values  using  the  above  procedure  is  given  as  a  function 


of  Rus  in  Fig.  7:  the  points  denote  the  experimental  values  and  the  curve 


is  obtained  from  the  digital  simulations.  The  observed  reasonable 


agreement  between  the  experimental  points  and  the  simulated  curve  supports 


the  underlying  validity  of  the  present  analysis. 


Consequently,  each  of  the  measured  k  .  (app)  values  in  Fig.  7  can 


be  used  to  estimate  k8,  (true)  for  this  system  with  reasonable  accuracy, 


by  evaluating  Rug  and  using  the  known  experimental  parameters 


(Cal,  P.  etc.)  to  simulate  a  series  of  corresponding  values  of  k  (app) 


and  kS,  (true),  as  exemplified  in  Fig.  4.  The  correct  value  of  k8  (true) 


will  be  that  corresponding  to  the  k8.  (app)  value  that  best  matches  the 


experimental  quantity.  Besides  providing  a  quantitatively  reliable 


estimate  of  k®b(true),  this  procedure  has  the  critical  virtue  of  yielding 
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an  assessment  of  the  maximum  rate  constant  that  can  be  evaluated  with  a 
given  instrument  and  experimental  technique. 

An  alternative  scheme  for  accounting  for  the  effects  of  solution 
resistance  and  amplifier  damping  would  be  to  correct  the  measured  AC 
currents  themselves.  However,  since  these  combined  effects  influence  the 
currents  in  a  complex  manner,  the  approach  followed  here  has  the  advantage 
of  greater  simplicity. 


Both  the  uncompensated  solution  resistance  and,  at  a  DME,  the 
amplifier  low-pass  filter  act  to  distort  AC  impedance  data  in  such  a  way 
to  make  the  apparent  observed  rate  constant  uniformly  smaller  than  the 
actual  value.  The  effects  of  these  two  factors  are  interactive  rather  than 
simply  additive,  and  always  combine  together  to  yield  an  error  in  k^b(app) 
that  is  larger  than  that  resulting  from  either  of  the  individual  factors 
alone.  Since  the  magnitude  of  the  distorting  effects  of  the  amplifier 
damping  depend  upon  the  total  solution  resistance  rather  than  the  uncom¬ 
pensated  portion  alone,  these  effects  will  be  retained  even  if  -  0. 

This  situation  might  be  regarded  as  presenting  a  bleak  picture  for  the 
evaluation  of  fast  electrode  kinetics,  at  least  using  conventional  AC 
impedance  measurements.  However,  the  same  simulations  that  are  used  to 
diagnose  the  presence  of  such  deleterious  effects  can  also  be  utilized  to 
improve  the  accuracy  by  which  moderately  fast  rate  constants  can  be 
evaluated . 


This  work  is  supported  by  the  Office  of  Naval  Research. 
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Equivalent  circuit  of  electrochemical  cell  employed  for  simulations. 

Fig,.  2 

S  8 

Logarithm  of  limiting  value  of  kQ^(app)  [i.e.  value  of  kQ^(app)  when 
s  s 

kQ^(true)  «  <*»]  ,  log  ko^(app,lim)  plotted  against  logarithm  of  uncompensated 
resistance,  log  R^g.  Simulation  conditions  are  C^i  -  5  pF  cm"^  and  p  - 
10  Q  cm.  Curve  A  refers  to  no  low-pass  filtering;  curves  B-D  for  filter 
time  constants  equal  to  0.03  s,  0.1  s,  and  for  a  pair  of  0.1  s  filters  in 
series . 

Fie.  3 

As  for  Fig.  2,  but  for  Cjji  -  20  pF  cm' ^ . 

FI  jy .  4 

g 

Plots  of  logarithm  of  apparent  (measured)  rate  constant,  log  k^Capp), 
against  logarithm  of  corresponding  "true"  value,  log  k*^(true),  for  Cdi  - 
20  pF  cm'2  and  p  -  500  0  cm.  Damping  conditions  for  curves  A-D  as  in  Fig. 
2. 

LU.  5 

The  cotangent  of  the  phase  angle  attributed  to  faradaic  current,  cot^, 
(see  text)  plotted  against  square  root  of  AC  frequency,  w1*,  for  various  R^ 
values,  with  C<n  -  20  p F  cm'^.  Straight  lines  shown  are  least  squares  best 
fits  to  points.  Key  to  R^  values:  o  ,  0  0;  v  ,  0.5  0;  O  ,  1  0;  A,  2  0. 
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Plots  of  log  kQ^(app)  versus  log  for  kQ^(true)  -  1  cm  s'1,  - 

20  pF  cm'^;  p  -  500  0  cm.  Damping  conditions  for  curves  A-D  as  in  Fig.  2. 

Fie.  7 

Experimental  values  of  k^ ( app )  for  Cp2Co+/°  in  dimethylsulfoxide 
(containing  0.1  K  tetrabutylammonium  hexafluorophosphate)  obtained  at  a  DME 
(solid  points)  as  a  function  of  the  uncompensated  solution  resistance, 
using  a  lock-in  amplifier  time  constant  of  0.1  s.  Solid  curve  is 
corresponding  simulated  data  obtained  for  k®^  (true)  -  0.7  cm  s'^.  (See 
text  for  further  details . ) 
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